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and non-Thermal Leptogenesis 
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Abstract: We consider a phenomenological extension of the minimal supersymmetric standard model 
(MSSM) which incorporates non-minimal chaotic inflation, driven by a quadratic potential in conjunction with 
a linear term in the frame function. Inflation is followed by a Peccei-Quinn phase transition, based on renormal- 
izable supeipotential terms, which resolves the strong CP and n problems of MSSM and provide masses lower 
than about 10 12 GeV for the right-handed (RH) (s)neutrinos. Baryogenesis occurs via non-thermal leptogenesis, 
realized by the out-of-equilibrium decay of the RH sneutrinos, which are produced by the inflaton's decay. 
Confronting our scenario with the current observational data on the inflationary observables, the light neutrino 
masses, the baryon asymmetry of the universe and the gravitino limit on the reheat temperature, we constrain 
the strength of the gravitational coupling to rather large values (~ 45 — 2950) and the Dirac neutrino masses 
to values between about 1 and 10 GeV. 
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L Introduction 

There is recently a wave of interest in implementing 
non-minimal chaotic inflation (nMCI) within both a non- 
supersymmetric (SUSY) [1-8] and a SUSY [9-13] frame- 
work. The main idea is to introduce a large non-minimal 
coupling of the infiaton field to the curvature scalar, 1Z. Af- 
ter that, one can make a transformation - from the Jordan 
frame (JF) to the Einstein (EF) one - which flattens the poten- 
tial sufficiently to support nMCI. The implementation of this 
mechanism within supergravity (SUGRA) has been greatly fa- 
cilitated after the developed [10] supercorformal approach to 
SUGRA. In particular, it is shown that the frame function can 
be related to a logarithmic type Kahler potential which en- 
sures canonical kinetic terms for the scalars of the theory and 
incorporates an holomorphic function, F, which expresses the 
non-minimal coupling of the infiaton field to 7Z. Until now, 
the proposed models [10-13] of nMCI within SUGRA are 
constructed coupling quadratically the infiaton superfield with 
another one in the superpotential - leading thereby to a quartic 
potential - and adopting a quadratic term for it in F. 

In this paper we propose a novel realization of nMCI within 
SUGRA, according to which the infiaton superfield is cou- 
pled linearly to another superfield in the superpotential of the 
model. As a consequence, a quadratic potential for the infiaton 
arises which supports nMCI, if the infiaton develops a linear 
coupling to 1Z. Actually, this set-up represents the SUSY im- 
plementation of the model of nMCI with n = — 1 introduced 
in Ref. [7]. In contrast to earlier models [10, 14] which re- 
lied on the same superpotential term - see also Ref. [15] -, no 
extra shift symmetry is imposed on the Kalher potential. The 
resulting mass of the infiaton lies at the intermediate scale and 
the inflationary observables are principally similar to those of 
nMCI with quartic - not quadratic - potential and therefore, in 
excellent agreement with the current observational data [16]. 

The inflationary model can be nicely embedded in a mod- 
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est phenomenological extension of the minimal supersymmet- 
ric standard model (MSSM) which incorporates a resolution 
of the strong CP problem [17] via a Peccei-Quinn (PQ) sym- 
metry. Note that there is an increasing interest [18, 19] in 
such models at present, since they provide us with two ad- 
ditional cold dark matter (CDM) candidates (axino and ax- 
ion) beyond the lightest neutralino. In our model, a PQ phase 
transition (PQPT), tied on renormalizable [11, 20] superpon- 
tential terms, can follow nMCI generating in addition, the fi 
term of MSSM and intermediate masses for the right handed 
(RH) [s]neutrinos, v% [Pf]. As a consequence, the light neu- 
trino masses can be explained through the well-known see- 
saw mechanism [21] provided that no large hierarchies occur 
in the Dirac neutrino masses. The possible formation [22] 
of disastrous domain walls can be avoided [23, 24] by intro- 
ducing extra matter superfields without jeopardizing the gauge 
unification of MSSM. The appearance of a Lagrangian quatric 
coupling of the infiaton ensures its decay to uf , whose the sub- 
sequent out-of-equilibrium decays can generate the Baryon 
Asymmetry of the Universe (BAU) via non-thermal leptoge- 
nesis (nTL) [25], consistent with the present data on neutrino 
data [26, 27]. Our model favors mostly quasi-degenerate v£ 
- as in Ref. [28] - which enhances the contribution from 
the self-energy corrections to leptonic asymmetries, without 
jeopardizing the validity of the relevant perturbative results, 
though. The constraints arising from BAU and the gravitino 
(G) limit [29-31] on the reheat temperature can be met pro- 
vided that the masses of G lie in the multi-TeV region. 

In Sec. II we present the basic ingredients of our model, 
Sec. Ill describes the inflationary scenario, and we outline the 
mechanism of nTL in Sec. IV. We then restrict the model pa- 
rameters in Sec. V and summarize our conclusions in Sec. VI. 
Throughout the text, the subscript of type , x denotes deriva- 
tion with respect to (w.r.t) the field \ (e.g., xx = d 2 /d\ 2 ); 
charge conjugation is denoted by a star and brackets are, also, 
used by applying disjunctive correspondence. 
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II. Model Description 

We focus on a PQ invariant extension of MSSM, which 
is augmented with (i) two superfields (P and P) which are 
necessary for the implementation of nMCI; (ii) three super- 
fields (S, $ and $) involved in the spontaneous breaking of 
the PQ symmetry, Z7(1)pq (iii) three RH neutrinos, v£, which 
are necessitated for the realization of the see-saw mecha- 
nism; (iv) n - to be determined below - pairs of SU (3)c 
triplets and antitriplets superfields, D a and D a respectively, 
(a = 1, n) in order to avoid the formation of domain walls 
- c.f. Ref. [23, 24] - and (v) an equal number of pairs of 
5[/(2)l doublet superfields, h a and h a in order to restore 
gauge coupling unification at one loop - see below. Be- 
sides the superfields in the points (iv) and (v), all the oth- 
ers are singlets under the Standard Model (SM) gauge group 
Gsm = SU(3) C x SU(2) L x U(1) Y . Besides the (color) 
anomalous [7(1)pq, the model also possesses an anomalous 
R symmetry U(1)r the baryon number symmetry U(1)b and 
two accidental symmetries U(1)d and U(l)h- The represen- 
tations under Gsm, and the charges under the global symme- 
tries of the various matter and Higgs superfields are listed in 
Table I. Note that the lepton number is not conserved in our 
model. 

In particular, the superpotential, W, of our model can be 
split into four parts: 



W = Wmssm + Wtw + Wcpq + W N r, 



(1) 



which are analyzed in the following: 

1 . Wmssm is the part of W which contains the usual 
terms - except for the [i term - of MSSM, supplemented by 
Yukawa interactions among the left-handed leptons and vf: 



Wmssm 



hr>ijd1QjH d 



huij u iQjHu 
<i Nij i 



:LjH u 



(2) 



Here, the group indices have been suppressed and summation 
over the generation indices i and j is assumed; the i-th genera- 
tion SU (2)l doublet left-handed quark and lepton superfields 
are denoted by Qi and Li respectively, and the SU{2)\ J sin- 
glet antiquark [antilepton] superfields by and [el and 
vf] respectively. The electroweak SC/(2)l doublet Higgs su- 
perfield, which couples to the up [down] quark superfields, is 
denoted by H u [Hd\. 

2. Wdw is the part of W which gives intermediate 
scale masses via ($) - see below - to D a — D a and h a — h a . 
Namely, 



Wdw = X D ^D a D a + \ ha <$>h a h a 



(3) 



Here, we chose a basis in the D a — D a and h a — h a space 
where the coupling constant matrices A^ a and A^ a are diag- 
onal. Although these matter fields acquire intermediate scale 
masses after the PQ breaking, the unification of the MSSM 
gauge coupling constants is not disrupted at one loop. In fact, 
if we estimate the contribution of D a , D a , and h a and h a to 
the coefficients b\, b 2 , and fe 3 , controlling [32] the one loop 
evolution of the three gauge coupling constants g\, g 2 , and g 3 , 
we find that the quantities b 2 — &i and 63 — b 2 (which are [32] 
crucial for the unification of g\,g 2 , and g 3 ) remain unaltered. 



TABLE I: Superfield Content of the Model 



Super- 
fields 



Representations 
under Gsm 



R 



Global Symmetries 
PQ B D 



Matter Fields 



Li 


(1,2,-1/2) 





-1 











et 


(1.1,1) 


2 


-1 













(1,1,0) 


2 


-1 











Qi 


(3,2,1/6) 


1 


-1 


1/3 








ut 


(3,1,-2/3) 


1 


-1 


-1/3 








dt 


(3,1,1/3) 


1 


-1 


-1/3 








Extra Matter Fields 


Da 


(3,1,-1/3) 


1 


1 





1 





D a 


(3,1,1/3) 


1 


1 





-1 





h a 


(1,2,1/2) 


1 


1 








1 


hu 


(1,2,-1/2) 


1 


1 








-1 




Higg 


s Fields 








H d 


(1,2,-1/2) 


2 


2 











H u 


(1,2,1/2) 


2 


2 











S 


(1,1,0) 


4 














$ 


(1,1,0) 





2 











$ 


(1,1,0) 





-2 











P 


(1,1,0) 


6 


1 











P 


(1,1,0) 


-2 


-1 












3. Wcpq is the part of W which is relevant for nMCI, 
the spontaneous breaking of f/(l) P Q, the decay of the infiaton 
and the generation of the masses of v ?'s and the /j, term of 
MSSM. It takes the form 



CPQ 



mPP + \ a S ($$ - Mpq) + \ iv c$ 



(4) 



where Mpq = / a /2with/ Q ~ (l0 10 - 10 12 ) GeV being the 
axion decay constant which coincides with the PQ breaking 
scale. The parameters A a and f a can be made positive by field 
redefinitions. From the terms in the right hand side (RHS) of 
Eq. (4) we note that the imposed symmetries disallow renor- 
malizable terms mixing P with some other superfields, which 
avoids undesirable instabilities faced in Ref. [11]. 

4. Wnr is the part of W which contains its non- 
renormalizable terms. Namely, we have 

Wm = K i^i + , P f±ii +Xr ^i, (5) 



nip 



TOp 



TOp 



where mp ~ 2.44 • 10 18 GcV is the reduced Planck scale. 
The first term in the RHS of Eq. (5) helps accomplish suffi- 
ciently low reheat temperature and leads to the production of 
uf's as dictated by nTL - see Sec. IV A. Finally, the third term 
provides the fi term of MSSM - see below. 

To get an impression for the role that each term in the RHS 
of Eqs. (3), (4) and (5) play, we display the SUSY potential, 
Ususy, induced from the following part of W 



Wei = Wcpq + W D w, 



(6) 
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which turns out to be 

Vsusy = m 2 (\P\ 2 + \P\ 2 ) + \\ iv ci>f + A a 5$| 2 
+ 4A^|^| 2 + A'|^-M| Q | 2 

+ | \ a S® + X Da D a D a + \ ha h a h a | 2 

+ \ 2 Da (\D a \ 2 + \D a \ 2 )\$\ 2 

+ AL (\h a \ 2 + \h a \ 2 ) |$| 2 , (7a) 

where the complex scalar components of the superfields 
P, P, S,<&,$, D a , D a ,h a , and h a are denoted by the same 
symbol as the corresponding superfields. From Eq. (6) and 
assuming [20] canonical Kahler potential for the hidden sec- 
tor fields, we can also derive the soft SUSY-breaking part of 
the inflationary potential which reads: 

V soit = m 2 Q <TC+ (rnBPP - a T A a SM 2 Q 
+ \ Da A Da $D a D a + \ ha A ha $h a h a 
+ \ a A a S<$>$ + \ iv cA iv c§vl + h.c.) (7b) 

where m^c , with 

Q = P, P, S, !>, $, vt, D ka , D ka , h la , h la (8) 

A a , Ai v c, Adh, Ah a , B and ax are soft SUSY-breaking mass 
parameters of order 1 TeV. From the potential in Eqs. (7a) 
and (7b), we find that the SUSY vacuum lies at 

(P) = (P) = = 0, (9a) 

(D ka ) = (D ka ) - (h la ) = (h la ) = 0, (9b) 

and 

(s) = |Aa ^ a |aTl , \(<h)\ = m)\ = m)\ = u, (90 

where the resulting (S) is of the order of TeV - cf. Ref. [20] 
- and we have introduced the canonically normalized scalar 
field 0$ = 2$ = 2$. Also, we use the subscripts k = 1, 2, 3 
and I = 1, 2 to denote the components of D a , D a and h a , h a , 
respectively. Note that, since the sum of the arguments of 
($), ($) must be 0, $ and $ can be brought to the real axis 
by an appropriate PQ transformation. After the spontaneous 
breaking of [7(1)pq, the third term in Eq. (4) generates in- 
termediate scale masses, for the vf's and, thus, seesaw 
masses [21] for the light neutrinos - see Sec. IV. The third 
term in the RHS of Eq. (5) leads to the fi term of MSSM, 

l — 1 2 

with ~ A^ ($) /mp, which is of the right magnitude if 
|($)| = f a /2 ~ 5 • 10 11 GeV, A M - (0.001 - 0.01). Fi- 
nally, since ($$) /m P = Mp Q /m P < m ~ 10 16 GeV - see 
Sec. IIIB and VB 1 - the second term in the RHS of Eq. (5) 
has no impact on our results. 

Nonetheless, Wei also gives rise to a stage of nMCI within 
SUGRA, if it is combined with a suitable Kahler potential, K, 
related to the frame function, Qci y i a 



In JF a specific form of flci's - see Ref. [10, 11] - en- 
sures canonical kinetic terms of the fields involved and a non- 
minimal coupling of the infiaton to 1Z represented by an holo- 
morphic function F(P). Going from JF to EF, and expanding 
the EF potential, V, along a stable direction - usually with all 
the fields besides infiaton placed at the origin - V takes the 
simple form 

Ucio - Vstjsy/7 {<?? , (ID 

where a = V%\P\ and / can be found expanding Oci- Van- 
ishing of the non-infiaton fields ensures, also, the elimination 
of some extra kinetic terms for scalars from the auxiliary vec- 
tor fields - see Ref. [9-11]. 

Let us emphasize here that the coupling of P to P is crucial 
in order to obtain the simple form of Vcio in Eq. (11), since 
only terms including derivatives of Wc\ w.r.t P survive in the 
EF SUGRA potential - see Sec. Ill A. This fact ensures the 
appearance of just one dominant power of a in the numera- 
tor of the SUGRA scalar potential. Such a construction is not 
possible, e.g., for a superpotential term of the form mP 2 . Ap- 
plying the strategy, described above Eq. (11) in our case, we 
can observe that along the direction 

= P = S = $ = <5> = i> c l = 

D ka = D ka = h la = h la = 0, (12a) 

with 9 = argP, Vsusy in Eq. (7a) becomes 

Ususy = \m 2 o 2 + A 2 M| Q . (12b) 

Clearly, for a 3> f a , Vsusy tends to a quadratic potential 
which can be flattened, according to Eq. (11), if / is mainly 
proportional to a, i.e., if F is a linear function of P with a 
sizable coupling constant c-jz. Therefore, we are led to adopt 
the following frame function: 

n cl = -3+^-^\P\ 4 -(F(P) + F*(P*)), G3a) 

with 0"'s defined in Eq. (8) and the non-minimal gravitational 
coupling 

F = 3c n P/V2m P , (13b) 

which breaks explicitly the imposed R and PQ symmetries 
during nMCI. In Eq. (13a) the coefficients kp and c-jz, for sim- 
plicity, are taken real. We remark that we add the third term 
in the RHS of Eq. (13a) to cure the tachyonic mass problem 
encountered in similar models [9, 10] - see Sec. Ill A. 

For P -C mp, we can show - see Sec. IV A - that an insta- 
bility occurs in the PQ system which can drive a PQPT which 
leads to the v.e.vs in Eq. (9c). Also, at the SUSY vacuum 
the explicit breaking of U(1)r x ?7(1)pq through Eq. (13b) 
switches off - see Eq. (9a). A closer look, however, reveals 
that instanton and soft SUSY breaking effects explicitly break 
U(1)r x {7(1)pq to Z2 x Z 2 ( n -6)^ as can be deduced from 
the solutions of the system 



K = -3mp ln(-0 C i/3). 



(10) 



4r = (mod27r) and 2(n-6)p- 12r = (mod 2tt) , (14) 
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where r and p are the phases of a U(1)r and J7(1)pq rotation 
respectively. Here, we take into account that the R charge 
of W and, thus, of all the soft SUSY breaking term is 4 and 
that the sum of the R [PQ] charges of the SU(3)c triplets and 
antitriplets is —12 [2(n — 6)]. Note that no loop-induced PQ- 
violating term - as this appearing in the first paper of Ref . [13] 
- is detected in our case. It is then important to ensure that 
Z 2 x Z 2 ( n _6) is not spontaneously broken by ($) and ($), 
since otherwise cosmologically disastrous domain walls are 
produced [22] during PQPT. This goal can be accomplished 
by adjusting conveniently the number n of D a — D a and h a — 
h a - see Table I. Indeed, when n = 5 or 7 we obtain 2p = 
(mod 2tt) and therefore, Z 2 x Z 2 ( n _ 6 ) is not spontaneously 
broken by (<p) and ($). The residual unbroken Z 2 subgroup 
of C/(1)pq can be identified with the usual matter parity of 
MSSM - see Table I - which prevents the rapid proton decay 
and ensures the stability of the lightest SUSY particle (LSP). 



m. The Inflationary Epoch 

In Sec. Ill A we describe the salient features of our infla- 
tionary model and in Sec. Ill B we extract the inflationary ob- 
servables. 



A. Structure of the Inflationary Potential 

The EF F-term (tree level) SUGRA scalar potential, Vcio, 
of our model is obtained from Wei in Eq. (6) and K in 
Eqs. (10) and (13a) by applying [9] 

Vcio = e K '™l (W^ - 3^) , (15a) 



with 



and 



Q7 



Fq = Wci,^ +A>Wci/m£, 



(15b) 



(15c) 



where the </>"'s are given in Eq. (8). From the resulting V"cio> 
we can deduce that along the field directions in Eq. (12a), 



Vcio = 



m 2 m 2 P xl + 4\lM* Q /m 4 P m 2 m %xl 



2/2 



2/2 



, (16a) 



where x a = a /nip and, according to the general recipe [9- 
11], the function 



/ = 1 + c n x a - xl/6 



(16b) 



expresses the non-minimal coupling of a to 1Z in JF. From 
Eq. (16a), we can verify that for c-jz 1 and x„ <C V®, 
Vci develops a plateau since Mpq <C mp - see Sec. IIIB. 
Along the trajectory in Eq. (12a), we can estimate the constant 
potential energy density 



2™2 



Vcio = 



ma rn'm 



2/2 



zc n 



(17a) 



and the corresponding Hubble parameter 



1/2 



Til 



V3mp V6c n 



(17b) 



In order to check the stability of the direction in Eq. (12a) 
w.r.t the fluctuations of the various fields, we expand them in 
real and imaginary parts according to the prescription 



ere" 



and X 



Xi + iX2 



where 



and 



V2 V2 ' 

X = P, S, #, Df, D ka , D ka , h ia , h la 



X=p,s, <p, <p, Vi, D ka , D ka , hi a , h ia , 
respectively. Along the trajectory in Eq. (12a) we find 
/ \ 



(K a p) = diag 




,2 / /,< 



J 2 , 

" " ' / 

\ 7+1 On elements/ 

2 _ 



(18a) 



(18b) 



(18c) 



(19a) 



1 J_. (19b) 

2 x a 



Consequently, we can introduce the EF canonically normal- 
ized fields, a, 9 and £, as follows - cf. Ref. [9-12]: 



1 



2 \ 1 / 2 2 \ 



where the dot denotes derivation w.r.t the JF cosmic time, t 
and the hatted fields are defined as follows 



|=J, *Q — J(jQ and ^ ^ 



(21) 



Note that 9 ~ Ja9 since Ja ~ y'i^mp - see Eq. (19b). 
On the other hand, we can show that during a stage of slow- 
roll nMCI, x — xl \fl since the quantity //2/ 3 / 2 x, involved 
in relating \ to x, turns out to be negligibly small compared 
with x- Indeed, the £'s acquire effective masses m£ 3> Hqi 
- see below - and therefore enter a phase of oscillations about 
X = with reducing amplitude. Neglecting the oscillating 
part of the relevant solutions, we find 



X-Xo 



/fe- 2 */ 3 and j~-2 X oV7ff. 



ci%f 



-2N/3 



, (22) 



where xo represents the initial amplitude of the oscillations, 
?>/3-ffci and we assume x(t = 0) = 0. Taking into 



Vx ..-y 



account the approximate expressions for a and the slow-roll 
parameter e, which are displayed in Sec. Ill B, we find 



2J3/2 



X 



m x 



2 — 2«X- 



(23) 
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TABLE II: The mass spectrum of the model during nMCI 



Fields 



Eingestates 



Mass Squared 



Bosons 



1 real scalar 

2 real scalars 

2 real scalars 
6 real scalars 
4 real 
scalars 
6n real scalars 
6n real scalars 
4n real scalars 
4n real scalars 



Pi,P 2 



Sl, S2 

V\i , V2i 

4>i±4>i 

-Dlfca, £>2ka 
Dlha,, Dlka. 
folia, fo-2ia 
folia, fo2ia 



dzm 2 x a /fj 2 ~ 4i?ci 
m 2 (2kpCjixl + (— c TC x CT + 
(6fcp - l)4z 2 )/2J 2 / 2 ) If 
(2A a Mpq/mp + m 2 z 2 ) /3/ 2 
(2A a M PQ /m 2 + m 2 z 2 ) // 2 
(m 2 4/3±A 2 M 2 Q /) // 2 

(2A a M PQ /m 2 + m 2 x 2 ) // 2 
(2A a M PQ /m 2 +m 2 xl)/f 
(2A a Mpq/mp + m 2 4) // 2 
(2A a M PQ /m 2 + m 2 * 2 ) // 2 



Fermions 



2 Weyl spinors 



x/2 



m 2 (6 + a ; 2 ) 2 /36/ 2 



The masses that the various scalars acquire during nMCI 
are presented in Table II. To this end, we expand Vcio in 
Eq. (15a) to quadratic order in the fluctuations around the di- 
rection of Eq. (12a). As we observe from the relevant eigen- 
values of the mass-squared matrices, no instability - as the 
one found in Ref. [1 1] - arises in the spectrum. In particular, 
it is evident that kp > 1 assists us to achieve positivity of the 
mass-squared associated with the scalars p 12 , mi - in accor- 
dance with the results of Ref. [9, 10]. It is remarkable that 
mass-squared corresponding to Df, Dk a , Dk a , fo; a , hi a are in- 
dependent of the relevant superpotential couplings A^ c , A_o a 
and Xha,- We have also numerically verified that the various 
masses remain greater than Hq\ during the last 50 — 60 e- 
foldings of nMCI, and so any inflationary perturbations of the 
fields other than the inflaton are safely eliminated. 

In Table II we also present the masses squared of chiral 
fermions of the model along the direction of Eq. (12a). In- 
serting these masses into the well-known Coleman- Weinberg 
formula [33], we can find the one-loop radiative corrections, 
V IC , in our model which can be written as 



V„ 



1 



ml 



64-7T 2 V 6 A 2 PA 2 s A 2 

m 2 m l m l 

+ 6mi In % + 2ml In -± + 2mi_ In -±- 

to- m % \ 

+ 20nm~ln-Hr A - 4 fe 

D yY 2 



4m ? ± ln A 2 



(24) 



where A = top /or is the cutoff scale of the effective the- 
ory - see Sec. V A - and the involved above masses squared 
to 2 , ml, to 2 , to 2 , to 2 , m - and ml are equal to the ones 

listed in the third column of Table II from top to the bottom. 



Note that the masses squared of all the extra matter fields are 
equal to to-. Based on the one-loop corrected EF potential 



Vci = Vcio + Kc 



(25) 



we can proceed to the analysis of nMCI in EE employing the 
standard slow-roll approximation [34]. It can be shown [35] 
that the results calculated this way are the same as if we had 
calculated them using the non-minimally coupled scalar field 
in JF. As expected and verified numerically, V TC does not af- 
fect the inflationary dynamics and predictions, in the major 
part of the allowed parameter space - see Sec. V B 1 - since 
the inflationary path already possesses a slope at the classical 
level - see below. 



B. The Inflationary Observables 

According to our analysis above, the universe undergoes a 
period of slow-roll nMCI, which is determined by the condi- 
tion - see e.g. Ref. [34]: 



max{e(cr), \rj(a)\} < 1, 



where 



TOt 



r = 



Vci,; 
Vci 



2J 2 



and 



m 



2 VqImS 
P 



Vci 



TOf, 



VW 

Vbi . 



VW 

Vbi 



4mp 

34^' 



Vci,. J,< 



(26a) 



Vci 



~ -4to p /3ck(t. 



(26b) 



Here, we employ Eqs. (17a) and (19b) and the following ap- 
proximate relations: 



Vci.o 



TO 2 TOp 



and Vci, 



2to 2 



1Z <T 



(27) 



The numerical computation reveals that nMCI terminates due 
to the violation of the ? criterion at a — uu which is calculated 
to be 



e(ai) = 1 



CTf = 



2 TOp 

y/3 c n 



(28) 



We note, in passing, that for a > Of the evolution of a - or a 
via Eq. (19b) - is governed by the equation of motion 



3H C 



da 
dt 



-Vci,* 



2\/2m 



(29) 



\ffdt. Using 



where t is the EF cosmic time with dt - 
Eqs. (26a) and (29), we can derive Eq. (23). 

The number of e-foldings, N*, that the scale fc* = 
0.002/Mpc suffers during nMCI can be calculated through 
the relation 



1 



A- ^1 

da — 

(ft Vci,s 



1 



da J' 



Vci 



(30) 
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where <r* [ct*] is the value of a [d] when fc* crosses the infla- 
tionary horizon. Given that cjf <C cr*, we can write ct* as a 
function of TV* as follows 



4mp 



V* — C7f) =>■ er» 



TOp- 



(31) 



The power spectrum of the curvature perturbations 
generated by a at the pivot scale fc* is estimated as follows 



27rm| 



6e~(<7*) 67rmpC7j' 



(32) 



where Eq. (31) is employed to derive the last equality of the 
relation above. Since the scalars listed in Table II are massive 
enough during nMCI, An can be identified with its central 
observational value - see Sec. V - with almost constant AT*. 
The resulting relation reveals that m is to be proportional to 
c-ji- Indeed we find 

m = 6irm P c n A n /N* m = 4.1 • 10 13 ct?, GeV, (33) 

for iV* ~ 55. At the same pivot scale, we can also calculate 
the (scalar) spectral index, n s , its running, a s , and the scalar- 
to-tensor ratio, r, via the relations: 



1 - 6e* + 2??* ~ 1 - 2/N*, 



16e* 



(n s - l) 2 ) - 2£, 
12/^ 2 , 



(34a) 
-2/iV 2 , (34b) 
(34c) 



where £ = mpVci,sVbi,s?s/Vci = m P VZerj i(T /J + 2-qe 
and the variables with subscript * are evaluated at a = a*. 
Comparing the results of this section with the observationally 
favored values [16], we constrain the parameters of our model 
in Sec. V B 1 . 



IV. Non-Thermal Leptogenesis 

A complete SUSY inflationary scenario should specify the 
transition to the radiation dominated era and also explain the 
origin of the observed BAU consistently with the G constraint. 
These goals can be accomplished within our set-up, as we de- 
scribe in this section. Namely, the basic features of the post- 
inflationary evolution are exhibited in Sec. IV A and the topic 
of nTL in conjunction with the present neutrino data is ana- 
lyzed in Sec. IV B. 



A. The General Set-up 

When nMCI is over, the inflaton continues to roll down to- 
wards the SUSY vacuum, Eqs. (9a), (9b) and (9c). Note that 
when x a < \/3\ a MpQ/m, one scalar originating from the 
superfields $ and $ - see Table I - acquires a negative mass 
squared triggering thereby the PQPT As the inflaton contin- 
ues its rolling, there is a brief stage of tachyonic preheating 



[36] which does not lead to significant particle production 
[37]. Soon afterwards, it settles into a phase of damped os- 
cillations about the minimum of the Vcio- Since no gauge 
symmetry is broken during nMCI, no superheavy bosons are 
produced and therefore no particle production via the mecha- 
nism of instant preheating [38] occurs. Also, since the inflaton 
cannot decay via renormalizable interactions to SM particles, 
effects of narrow parametric resonance [36] are also absent in 
our regime. 

Nonetheless, the standard perturbative approach to the in- 
flaton decay provides a very efficient decay rate. Namely, at 
the SUSY vacuum the fields involved acquire the v.e.vs shown 
in Eqs. (9a), (9b) and (9c) giving rise to the mass spectrum 
presented in Table III. There we can show the mass, mi, of the 
(canonically normalized) inflaton P and the masses M^c of 
the RH [s]neutrinos, vl [Df], which play a crucial role in our 
scenario of nTL. Note that since ($) = ($) = Mpq <C mp, 
(f2) ~ —3 and so (/) ~ 1. Therefore, apart from P, the EF 
canonically normalized field are not distinguished from the 
JF ones at the SUSY vacuum. On the other hand, P can be 
expressed as a function of P through the relation 



P 
P 



(J) where (J) = Jl + 3c n /2. 



(35) 



Making use of Eq. (33) we can infer that mi is kept inde- 
pendent of c-ji and almost constant at the level of 10 13 GeV. 
Indeed, 



mi ~ — = d- — ~ 2V3 7 rm P ^ r 
(J) V 3c R N 



10 13 GeV, (36) 



where the WMAP7 value of A-jz - see Sec. V A - is employed 
in the last step of the relation above. In the expressions of the 
various eigenstates listed in Table III, we adopt the following 
abbreviations 



5$ = $ - Mp Q , 5$ = $ - M] 



and 



(37a) 



(37b) 



where ip x with x = P,P, S, $, $, D ka , D ka , hi a , and /i; a 
denote the chiral fermions associated with the superfields 
P, P, S, <&, $, D ka , D ka , hi a , and hi a respectively. The eigen- 
states ^_ and <5$_, with 



5$± = (<5$ ± <5<I>) /V2, 



(38) 



contain the components of the axion supermultiplet. Namely 
axion [saxion] can be identified with the phase [modulus] of 
the complex field whereas ip_ can be interpreted as the 
axino. Note that the zero masses of saxion and axino can 
be replaced with masses of order 1 TeV if we take into ac- 
count the soft SUSY breaking masses - see discussion below 
Eq. (9c). 

The decay of P commences when mj becomes larger than 
the expansion rate and is processed via the first coupling in 
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TABLE III: The mass spectrum of the model at the SUSY vacuum 



Eigenstates 


Eigenvalues 


Scalars 


Fermions 


(Masses) 


P,P 


(V^p ± Vp) /y/2 


mi = m/(J) 


S, (<5$ + <5$) /y/2 


(i>s ± V+) /V2 


mpQ = a/2A M P q 


(<5$ - <5<E>) /y/2 


(V* - V*)/\/2 





Of 




Mi^c = 2Ai„cM PQ 


Dka,, D k& 


^*.> V>S fca 


m. Da = A_o a MpQ 


hla, hla 


V^a> V>h ia 


mf. a = A fca M P q 



the RHS of Eq. (5), into S and Of and 5, i>? and S<S>+ or <5$_. 
The relevant Lagrangian sector is 

C dc = -^X^Svf ( Mpq + S *+- S *A +h . c . (39) 

which arises from the cross term of the F-term, corresponding 
to P, of the SUSY potential derived from the superpotential 
terms in Eqs. (4) and (5). Note that we have no c^-induced 
decay channels as in Ref. [12], since (P) = 0. The interaction 
above gives rise to the following decay width 




where we take into account that mi ^> mpQ and mi Mj U <= . 
These prerequisites are safely fulfilled when A a and A^c re- 
main perturbative, i.e. A a , A^c < - see Table III. From 
the two contributions to Ti, the dominant one is the second 
one - the 3 -body decay channel - originating from the two 
last terms of Eq. (39). 

Taking also into account that the decay width of the pro- 
duced Of, Y iv c , is much larger than Tj- see below - we can in- 
fer that the reheat temperature, T r h, is exclusively determined 
by the P decay and is given by [39] 

T rll =(^-) 1/4 VTW, (41) 

where g„ ~ 232.5 counts the effective number of the relativis- 
tic degrees of freedom at temperature T r h for the (s)particle 
spectrum of MSSM plus the particle content of the axion su- 
permultiplet. Although the factor before the square root of 
Eq. (41) differs [39] slightly from other calculations of T r h - 
cf. Ref. [25] - the numerical result remains pretty stable and 
close to 10 8 GeV - see Sec. V B 1 . 

If T r h <C Mi„ c , the out-of-equilibrium condition [40] for 
the implementation of nTL is automatically satisfied. Subse- 
quently, Of decay into H u and Li or H* and L* via the tree- 
level couplings derived from the second term of the second 
line of Eq. (2). Interference between tree-level and one-loop 
diagrams generates a lepton-number asymmetry (per vf de- 
cay) £j [40], when CP is not conserved in the Yukawa cou- 
pling constants h^ij - see Eq. (2). The resulting lepton- 
number asymmetry after reheating can be partially converted 



through sphaleron effects into baryon-number asymmetry. 
However, the required T rh must be compatible with con- 
straints for the G abundance, Yg, at the onset of nucleosyn- 
thesis (BBN). In particular, the B yield can be computed as 

Y B = -0.35^— V Br i£j with Br, = =^-= (42) 



the branching ratio of P to Of - see Eq. (40). In the for- 
mula above the first numerical factor (0.35) comes from the 
sphaleron effects, whereas the second one (5/4) is due to the 
slightly different calculation [39] of T rh - cf. Ref. [25]. On 
the other hand, the G yield due to thermal production at the 
onset of BBN is estimated to be [3 1 ] : 

F g ~ 1.9 • 10" 22 T rh /GeV. (43) 

where we assume that G is much heavier than the gauginos. 
Let us note that non-thermal G production within SUGRA 
is unlikely in our scenario, since these contributions are [41] 
usually proportional to the v.e.v of the inflaton which is zero 
in our case. 

Both Eqs. (42) and (43) calculate the correct values of the 
B and G abundances provided that no entropy production oc- 
curs for T < T r h - see also Sec. V A. This fact can be eas- 
ily achieved within our setting. Indeed, following the argu- 
ments of Ref. [11], one can show that the PQ system com- 
prised of the fields S and decays via the third term in 
the RHS of Eq. (5) before its domination over radiation, for 
all relevant values of Aj's. Regarding the saxion, we 
can assume that it has mass of the order of 1 TeV, its decay 
mode to axions is suppressed (w.r.t the ones to gluons, hig- 
gses and higgsinos [19, 42, 43]) and the initial amplitude of 
its oscillations is equal to f a . Under these circumstances, it 
can [42] decay before domination too, and evades [43] the 
constraints from the effective number of neutrinos for the / a 's 
and T r h's encountered in our model. As a consequence of its 
relatively large decay temperature, the LSPs produced by the 
saxion decay are likely to be thermalized and therefore, no 
upper bound on the saxion abundance is [43] to be imposed. 
Finally, if axino is sufficiently light it can act as a CDM can- 
didate [18, 19] with relic abundance produced predominantly 
thermally - due to the relatively large T r h- Otherwise, it may 
enhance [19] non-thermally the abundance of a higgsino-like 
neutralino-LSP, rendering it a successful CDM candidate. 



B. Lepton-Number Asymmetry and Neutrino Masses 

As mentioned above, the decay of Of, emerging from the 
P decay, can generate a lepton asymmetry, e%, caused by the 
interference between the tree and one-loop decay diagrams, 
provided that a CP- violation occurs in h^i/s. The produced 
£i can be expressed in terms of the Dirac mass matrix of Vi, 
771d, defined in a basis (called vf -basis henceforth) where vf 



V Constraining the Model Parameters 



8 



are mass eigenstates, as follows: 



where 



Ei^ 7 lm Kwd)-, (xij,Vi, Vj) + F v (xij) 



8ir(J/„»J)mD)i 



where we take (H u ) ~ 174 GeV, for large tan (3 and 



M jv c r 

and y { = 



M,, 



(WpMp),, 
Sn(H u ) 2 



(44a) 



(44b) 



Also Fy and Fs represent, respectively, the contributions 
from vertex and self-energy diagrams which in SUSY theo- 
ries read [44-46] 



Fy(x) 



•lull 



(44c) 



and 



F s (x, y, z) = 



-2x{x 2 - 1) 



(x 2 — 1 — x 2 z In x 2 /ir) 2 + (x 2 z — y) 2 ' 



(44d) 

with the latter expression written as given in Ref. [46]. When 



A iH > 1 and Am > 1 with A 



ijk 



cl - II 



(45) 



(no summation is applied over the repeated indices) we can 
simplify Fs expanding it close to x ~ 1 as follows 



2x^1 1 
l-x 2 ~ 1-x ~ 2 



(46) 



The involved in Eq. (44a) mp can be diagonalized if we define 
a basis - called weak basis henceforth - in which the lepton 
Yukawa couplings and the SU{2)\, interactions are diagonal 
in the space of generations. In particular we have 

U^muU^ = d D = diag (mm, m 2 D, m 3D ) , (47) 

where U and U c are 3 x 3 unitary matrices which relate Li 
and v£ (in the uf -basis) with the ones L\ and vf in the weak 
basis as follows: 



L 1 = LU and v cl = U c v c . 



(48) 



Here, we write LH lepton superfields, i.e. SU(2)l doublet 
leptons, as row 3-vectors in family space and RH anti-lepton 
superfields, i.e. SU{2)\ J singlet anti-leptons, as column 3- 
vectors. Consequently, the combination m^mo appeared in 
Eq. (44a) turns out to be a function just of g?d and U c . Namely, 



m^TOD = U cf dld D U c . 



(49) 



The connection of the leptogenesis scenario with the low 
energy neutrino data can be achieved through the seesaw for- 
mula, which gives the light-neutrino mass matrix m u in terms 
of TOiD and Mi„c. Working in the v\ -basis, we have 



m v = -m D d v } ra D , 



(50) 



d v c =d\ag(M lv c,M 2 „c,M 3v c) 



(51) 



with Mi u c < M 2u c < M 3u c real and positive. Solving 
Eq. (47) w.r.t tod and inserting the resulting expression in 
Eq. (50) we extract the mass matrix 



fh v = U^m v U* = -d u U c d-}U cJ d u , 



(52) 



which can be diagonalized by the unitary PMNS matrix satis- 
fying 



fh v = U* d\ag(m lv ,m 2v ,m 3v ) W v 
and parameterized as follows: 

-is' 



(53) 



C12C13 S12C13 s^e 

U21p U 22v S23C13 
U 3 \ u U 32v C23C13 



Here 



V 2 \v = ~c 23 s 12 - s 23 c 12 s 13 e 



iS 



U 22v = c 23 c 12 - s 2 3Si 2 Si 3 e 

u 32v 



S23S12 - c 23 ci 2 s 13 e 



iS 
JS 



-S23C12 - c 23 si 2 si 3 e 



(54) 



(55a) 
(55b) 
(55c) 
(55d) 



with Cij := cos 6ij, := sin^ and S the CP-violating 
Dirac phase. The two CP-violating Majorana phases <pi and 
ip 2 are contained in the matrix 



V = diag (e- tVl / 2 ,e- l ^/ 2 ,lj . 



(56) 



Following a bottom-up approach, along the lines of 
Ref. [47], we can find fh v via Eq. (53) using as input pa- 
rameters the low energy neutrino observables, the CP violat- 
ing phases and adopting the normal or inverted hierarchical 
scheme of neutrino masses. Taking also m^o as input param- 
eters we can construct the complex symmetric matrix 



W = -d^-fh v d^ 1 = U c d u .U cJ (57) 

- see Eq. (52) - from which we can extract d u c as follows: 

d-J = U^WW^U . (58) 

Note that WW^ is a 3 x 3 complex, hermitian matrix and 
can be diagonalized following the algorithm described in 
Ref. [48]. Having determined the elements of U c and the 
M iu c's we can compute m D through Eq. (49) and the £j's 
through Eq. (44a). 

v. Constraining the Model Parameters 

We exhibit the constraints that we impose on our cosmolog- 
ical set-up in Sec. V A, and delineate the allowed parameter 
space of our model in Sec. V B. 
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A. Imposed Constraints 



differences: 



The parameters of our model can be restricted once we im- 
pose the following requirements: 

1. According to the inflationary paradigm, the hori- 
zon and flatness problems of the standard Big Bang cosmol- 
ogy can be successfully resolved provided that TV* defined by 
Eq. (30) takes a certain value, which depends on the details of 
the cosmological scenario. Employing standard methods [7], 
we can easily derive the required TV* for our model, consistent 
with the fact that the PQ oscillatory system remains subdomi- 
nant during the post-inflationary era. Namely we obtain 



~ 22.5 + 2 In 
3 1 GcV 



VciM 1/4 
1 GcV 

■ + ii n /M. 

2 /(a*) 



4 VcM) 1/4 
3 n lGcV 



(59) 



2. The inflationary observables derived in Sec. Ill B are 
to be consistent with the fitting [16] of the WMAP7, BAO 
and Ho data. As usual, we adopt the central value of An, 
whereas we allow the remaining quantities to vary within the 
95% confidence level (c.l.) ranges. Namely, 



A TC ~4.93-10" 5 , 
0.944 < n s < 0.992, 
-0.062 <a s < 0.018, 
r < 0.24. 



(60a) 
(60b) 
(60c) 
(60d) 



3. For the realization of nMCI, we assume that c-jz takes 
relatively large values - see e.g. Eq. (16a). This assumption 
may [5, 49] jeopardize the validity of the classical approxi- 
mation, on which the analysis of the inflationary behavior is 
based. To avoid this inconsistency - which is rather question- 
able [10, 49] though - we have to check the hierarchy between 
the ultraviolet cut-off scale [7], A = m-p/cn, of the effec- 
tive theory and the inflationary scale, which is represented by 
V& (v*) 1 ^ or > l ess restrictively, by the corresponding Hub- 
ble parameter, H* — Vci( cr *) 1 ^ 2 / V3mp. In particular, the 
validity of the effective theory implies [49] 



(a) T>ciM 1/4 < A or (b) £T* < A. 



(61) 



4. To ensure that the inflaton decay according to the 
lagrangian part of Eq. (39) is kinematically allowed we have 
to impose the constraint - see Table III: 



mi > 2m PQ + Mj, 



2m PQ + M iv c < 10 1J GeV, (62) 



where we make use of Eq. (36). This requirement can be 
easily satisfied by constraining A a and \ v c to values lower 
than the perturbative limit. As the inequality in Eq. (62) gets 
strengthened, the accuracy of Eq. (40) where masses of the 
decay products are neglected, increases. 

5. From the solar, atmospheric, accelerator and reactor 
neutrino experiments we take into account the following in- 
puts [26] - see also Ref. [27] - on the neutrino mass-squared 



Ami = (7.59t°'? s ) • 10" 



eV 2 



'21 — V'- UJ -0.18; 
Ami = (2.5±g : ?g[-2.4±g;8a) 

on the mixing angles: 



0.016 



sin 2 


#12 


- 31 2+ 0017 


sin 2 


#13 


- 013+ 007 


sin 2 


#23 


- 52+ - 06 



10" d eV" 



+0.008] 
-0.006J ) 



[0.52 ±0.06], 
and on the CP-violating Dirac phase: 

* = - (o.6il8:3 [o.4i±r 7 5 ])- 



(63a) 
(63 b) 

(63c) 
(63d) 
(63e) 

(63f) 



for normal [inverted] neutrino mass hierarchy. In particular, 
rriivS can be determined via the relations 



m 2v = \J m\ v + Am. 



and 



m3v — \/ "ii/ 
for normally ordered (NO) m^'s or 

m lv = \jm\ v + | Am 2 . 



31 



311 



(64a) 



(64b) 



(64c) 



for invertedly ordered (IO) m^'s. The sum of m^'s can be 
bounded from above by the WMAP7 data [16] 



< 0.58 eV 



(65) 



at 95% c.l. This is more restrictive than the 95% c.l. upper 
bound arising from the effective electron neutrino mass in /?- 
decay [50]: 



m fi ■= \EiULm iv \ < 2.3 eV. 



(66) 



However, in the future, the KATRIN experiment [51] expects 
to reach the sensitivity of mp ~ 0.2 eV at 90% c.l. 

6. The interpretation of B AU through nTL dictates [ 1 6] 
at 95% c.l. 



Y B = (8.74 ±0.42) • 10" 



(67) 



7. In order to avoid spoiling the success of the BBN, an 
upper bound on Yg is to be imposed depending on the G mass, 
mg, and the dominant G decay mode. For the conservative 

case where G decays with a tiny hadronic branching ratio, we 
have [31] 



io- 14 

2.5- 10~ 14 
4.3- 10~ 14 
HT 13 



for nig 



0.69 TeV 
5TcV 
8 TeV 
10.6 TeV. 



(68) 



As we see below, this bound is achievable within our model 
model only for > 8 TeV. The bound above may be some- 
how relaxed in the case of a stable G. 
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1 10 1 10 

m(10 16 GeV) m(10 16 GeV) 

(a) (b) 

FIG. 1: The allowed by Eqs. (59), (60a), (600) and (61b) values of c-jz (solid line), mi - given by Eq. (36) ~ (dashed line) and T rb - given 
by Eq. (41) - (dot-dashed line) [ot (solid line) and er» (dashed line)] versus m (a) [(b)] for Ai <C A2 = A3 = 0.5. The light gray and gray 
segments denote values of the various quantities satisfying Eq. (61a) too, whereas along the light gray segments we obtain ct» > mp . Values 
of the parameters to the right of the lined region correspond to n B 's lying within its 68% c.l. observationally favored region. 



B. Results 

As can be easily seen from the relevant expressions in 
Sees. II and IV B, our cosmological set-up depends on the fol- 
lowing independent parameters: 

to, A a , A M , kp, A^ f a , n, mi v , m jD , tpi and</? 2 , 

where mt> v is the low scale mass of the lightest of iVs and 
can be identified with m\ v [m 3l/ ] for NO [IO] neutrino mass 
spectrum. We do not consider c-ji and Xi U a as independent pa- 
rameters since cu is related to to via Eq. (33) while X iu a can 
be derived from the last six parameters above which affect ex- 
clusively the Yl calculation and can be constrained through 
the requirements 5 and 6 of Sec. V A. Note that the X iu o 's can 
be replaced by M iu a 's given in Table III keeping in mind that 
perturbativity requires X iv a < \f^K<yrM iv c < 3.5/ a . Recall 
also that V TC in Eq. (24) is independent of Xd&, ^ha S> A a 
and depends only on n, which is set equal to 5 for definite- 
ness. To facilitate the realization of see-saw mechanism, we 
take f a = 10 12 GeV. This choice makes also possible the 
generation of the /i term of MSSM through the PQ symme- 
try breaking, since /r ~ 1 TeV is obtained for A M = 0.01, 
whereas lower / a 's dictate larger A M 's. Moreover, our com- 
putation reveals that e\ in Eq. (44a) is mostly smaller than 
£2 and £3. Therefore, fulfilling the baryogenesis criterion en- 
forces us to consider Bri <C Br 2 ,3 or Ai <C A 2 ,3 ~ 0.1. Since 
e 2 and £ 3 are of the same order of magnitude, the resulting 
Yb does not depends crucially on A2/A3. Therefore we be- 
lieve that A 2 = A 3 = 0.5 is a representative choice - e.g., 
we explicitly checked that the option A 2 = 0.1 and A 3 = 0.9 
or A2 = 0.9 and A3 = 0.1 lead to similar results. Finally, 
our results are independent of A a and kp provided Eq. (62) is 
fulfilled and the positivity of to| - see Table II - is ensured, 
respectively. To facilitate the achievement of these objective, 
we get A a = 0.01 and kp = 1. 



Summarizing, we set throughout our calculation: 

kp = 1, Ai < 0.01, A 2 = A 3 = 0.5, n = 5, (69a) 
A M = A a = 0.01 and f a = 10 12 GeV. (69b) 

The selected values for the above quantities give us a wide and 
natural allowed region for the remaining fundamental parame- 
ters of our model, as we show below concentrating separately 
in the inflationary period (Sec. V B 1 ) and in the stage of nTL 
(Sec. VB2). 



1. The Stage of non-Minimal Inflation 

For nMCI, we use as input parameters in our numerical 
code o~ * , to and Cfc. For every chosen c-ji > 1 we restrict m 
and cr* so that the conditions in Eq. (59) - with T r h evaluated 
consistently using Eq. (41) - and (60a) are satisfied. Let us 
remark that, in our numerical calculations, we use the com- 
plete formulae for Vqi - see Eq. (25) -, N*, the slow-roll 
parameters and Ar in Eqs. (30), (26a), (26b), (32) and not 
the approximate relations listed in Sec. Ill B for the sake of 
presentation. 

Our results are displayed in Fig. 1, where we draw the al- 
lowed values of c-n (solid line) m\ (dashed line)and T r h (dot- 
dashed line) [erf (solid line) and cr* (dashed line)] versus m 
- see Fig. l-(a) [Fig. l-(b)]. The constraint of Eq. (61b) is 
satisfied along the various curves whereas Eq. (61a) is valid 
only along the gray and light gray segments of these. Along 
the light gray segments, though, we obtain 0% > mp. The 
lower bound on to is derived from the saturation of the up- 
per bound of inequality in Eq. (60b) whereas the upper bound 
comes from the fact that the enhanced resulting to's destabi- 
lize the inflationary path through the radiative corrections in 
Eq. (25) - see Eq. (24). Indeed, V YC starts to influence the 
inflationary dynamics for m > 1.5 • 10 16 GeV, and conse- 
quently, the variation of cjf as a function of c-ji or m - drawn in 
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Fig. l-(b) - deviates from the behavior described in Eq. (28). 
On the contrary the variations of cr* follows Eq. (31). 
In all, we obtain 

Tfl 

45<c^< 2950 and 2.5 < < 102 (70) 

10 10 GeV 

for N* ~ 54.5. From Fig. l-(a), we observe that m depends 
on c-jz almost linearly whereas mi remains close to 10 13 GeV 
as we anticipated in Eqs. (33) and (36), respectively. As a 
result of the latter effect, T r h given by Eq. (41) remains also 
almost constant. As m (or c-jz) decreases below its maximal 
value in its allowed region in Eq. (70), we obtain 

0-965 < n s < 0.991, (71a) 
6.5 < -a s /10- 4 < 12, (71b) 
3.1 < r/10" 3 < 7.3. (71c) 

Clearly, the predicted n s , a s and r can lie within the allowed 
ranges given in Eqs. (60b), (60c) and (60d) respectively. In 
particular, values of the various parameters plotted in Fig. 1, 
which lie to the right of the lined regions correspond to n s ~ 
(0.965 - 0.98). This result is consistent with the 68% c.l. 
observationally favored region - see Eq. (60b). It is notable, 
however, that n s increases impressively for <r*/mp > \/6, 
contrary to the situation in models of nMCI with quadratic 
coupling to 1Z where n s remains constantly close to its central 
observational favored value in Eq. (60b) - cf. Ref. [11]. 
As regards the G abundance, employing Eq. (43), we find 

3.5 < Yg/10- 14 < 8.4 (72) 

as m varies within its allowed range in Eq. (70). Compar- 
ing this result with the limits of Eq. (68), we infer that our 
model is consistent with the relevant restriction for rag > 
(8 - 10) TeV. 



TABLE IV: Parameters yielding the correct BAU for various neu- 
trino mass schemes. 



Parameters 


Cases 




A 


B 


C 


D 


E 


F 


G 




Normal 
Hierarchy 


Degenerate 
Masses 


Inverted 
Hierarchy 



Low Scale Parameters 



m lv /0.1 eV 


0.05 


0.1 


0.5 


1. 


0.7 


0.5 


0.49 


m 2 „/0.1 eV 


0.1 


0.13 


0.51 


1.0 


0.705 


0.51 


0.5 


m 3 „/0.1 eV 


0.5 


0.51 


0.71 


1.12 


0.5 


0.1 


0.05 


^m^/O.leV 


0.65 


0.74 


1.7 


3.1 


1.9 


1.1 


1 


mp/0.1 eV 


8 • 10" 3 


0.013 


0.19 


0.46 


0.3 


0.42 


0.44 


fi 


7T 


7T 





7r/4 


tt/4 


7r/4 


tt/4 


<P2 








5tt/6 


7T 


7T 


tt/2 


tt/4 


Leptogenesis-Scale Parameters 


miD/GeV 


2 


2.5 


4 


8 


9 


6 


5 


TO2D /GeV 


3 


3.49 


5 


9.3 


6 


3 


1 


m 3D /GeV 


6.7 


4 


8 


11 


4.7 


2 


2.1 


Mi^c/10 11 GeV 


2.5 


2.4 


3.3 


6.5 


4.6 


1 


0.3 


M 2 „c/10 n GeV 


11 


7.3 


5.2 


8.13 


4.9 


5.56 


4.3 


Ma^/lO 11 GeV 


17 


7.6 


6 


8.36 


8.6 


6.7 


5.1 


A^/10 4 


4.5 


0.5 


0.96 


0.05 


0.3 


1.6 


1.4 


AW10 4 


1.6 


0.54 


0.48 


0.04 


0.3 


1.2 


3.5 



(with i = 2 and j = 3 except for case E where i = 1 and j = 2) 



Resulting B- Yield 





8.3 


7.4 


6.3 


3.3 


7.2 


9.3 


4.8 


io n y s 


8.7 


8.85 


8.98 


8.4 


8.9 


8.96 


8.95 



2. The Stage of non-Thermal Leptogenesis 

As we show above, the stage of nMCI predicts almost con- 
stant values of mi and T r h - recall that we consider Aj's of 
the order of 0.1. In other words, the post-inflationary evolu- 
tion in our set-up is largely independent of the precise value 
of m in the range of Eq. (70). As a consequence, Yb cal- 
culated by Eq. (42) does not vary with m, contrary to the 
naive expectations. Just for definiteness we take through- 
out this section m = 4.2 • 10 15 GeV which corresponds 
to c n = 100, n s = 0.969, mi = 3.4 • 10 13 GeV and 
T rh = 2.1 • 10 8 GeV (Yq ~ 4 • 10~ 14 ) - recall that we use 
Ai < A 2 = A 3 = 0.5. 

On the contrary, Yb in our approach depends crucially on 
the low energy parameters related to the neutrino physics. In 
our numerical program, for a given neutrino mass scheme, 
we take as input parameters: m&,,<£i,v?2 and the best-fit 
values of the neutrino parameters listed in the paragraph 5 
of Sec. VA. We then find the renormalization group (RG) 
evolved values of these parameters at the scale of nTL, Ax,, 
which is taken to be Ax = mi, integrating numerically the 



complete expressions of the RG equations - given in Ref. [52] 
- for m iu , 8ij, 5, ipi and ip 2 - In doing this, we consider the 
MSSM with tan /3 ~ 50 (favored by the preliminary LHC 
results [53]) as an effective theory between Ax and a SUSY- 
breaking scale, Msusy = 1.5 TeV. Below Msusy the run- 
ning of the various parameters is realized considering the par- 
ticle content of SM with a mass of about 120 GeV for the 
light Higgs. Following the procedure described in Sec. IV B, 
we evaluate M iu c at Ax taking m^ as free parameters. In our 
approach we do not consider the running of m^ and M iv a 
and therefore we give their values at Ax- 

We start the exposition of our results arranging in Table IV 
some representative values of the parameters leading to the 
correct BAU for normally hierarchical (cases A and B), de- 
generate (cases C, D and E) and invertedly hierarchical (cases 
F and G) neutrino masses. For comparison we display the 
B-yield with (Yb) or without (Yg) taking into account the 
RG effects. We observe that the two results differ apprecia- 
bly especially in the cases with degenerate or IO m^'s. As 
it is evident from the m^o 's chosen, our model is not com- 
patible with any GUT-inspired pattern of large hierarchy be- 
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FIG. 2: Contours in the m,2r> — mso (a) Mi„c — M 2v c (b) and M 2 v" — M 3l/ c (c) plane yielding the central Yb in Eq. (67), for various 
(m tl/ , mm ,<pi,(p2)'s indicated next to the graph (a) and NO [IO] m iv 's (black [gray] lines). 



tween the TOjd's. In particular, we need mm < ™ 2 d < m 3 D 
< m2D < ™id] for NO [IO] m^'s (cases A, B, C and 
D [cases E, F and G]). 

From Table IV we also notice that the achievement of Yb 
within the range of Eq. (67) dictates mostly proximity be- 
tween two of the Mj„c's. Indeed, except for the case A, we 
obtain M 2v c/M lv c ~ 1.06 in case E and M Sv c/M 2v c < 1.2 
in the residual cases. However, it is clear from the displayed 
Ajjj's and A^j's (with i = 2 and j ' = 3 for all the cases be- 
sides case E where i = 1 and j = 2) that in our framework 
the conditions of Eq. (45) are comfortably retained and there- 
fore, our proposal is crucially different from that of resonant 
leptogenesis [44—46] - it rather resembles that of Ref. [28] . 
On the other hand, the correctness of Yb in the case A entails 
M 2 „c and M 3l/ c [X 2 ^ and X^] roughly larger than 10 12 GeV 
[unity]. In all cases the current limit of Eq. (65) is safely met 
- the case D approaches it -, while mp turns out to be well 
below the projected sensitivity of KATRIN [51]. 

To highlight further our conclusions inferred from Table IV, 
we can fix (mi„ for NO m^'s or m^ v for IO m^'s) 
miD, <fi and ip 2 to their values shown in this table and vary 
to 2 d an d m 3D so that the central value of Eq. (67) is achieved. 
The resulting contours in the TO2D — plane are presented 
in Fig. 2-(a) - since the range of Eq. (67) is very narrow the 



possible variation of the drawn lines is negligible. The re- 
sulting values of Mj v c are displayed in M lu ^ — M 2v " and 
M 2v c — M^c plane - see Fig. 2-(b) and Fig. 2-(c) respec- 
tively. The conventions adopted for the types and the color of 
the various lines are also described next to the graphs (a) of 
Fig. 2. In particular, we use black [gray] lines for NO [IO] 
mi v 's. Besides the case with m\ v = 0.005 eV we observe 
that every curve in all graphs has two branches and not large 
hierarchies allowed in the sectors of both the tod's and M„c's. 
Note that the black contour for m\ v = 0.01 eV in Fig. 2-(c) 
is included within the one for m,\ v — 0.1 eV and so, it is not 
quite distinguishable. For m\ v — 0.005 eV, saturates the 
perturbation limit. Since we expect that the A,„c's increase 
[54] due to their RG running from low to higher scale, our 
results do not jeopardize the validity of the conventional per- 
turbation approach up to the scale A^. In all cases we find 
1 < m tD /GeV < 10. 

It is worth emphasizing that, although our mechanism of 
nTL is connected with the specific inflationary model under 
consideration, it can have a much wider applicability. It can be 
realized within other models of inflation with similar inflaton 
mass and reheat temperature, since it is largely independent of 
the details of the inflationary phase but restricts mainly the yet 
unknown parameters of neutrino physics (mi v ,miT>,(pi,(p 2 ). 
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vi. Conclusions 

We investigated a novel inflationary scenario in which the 
inflaton field appears in a bilinear superpotential term and in a 
linear holomorphic function included in a logarithmic Kahler 
potential. The latter function can be interpreted in JF as a 
non-minimal coupling to gravity, whose the strength is con- 
strained so as the EF inflationary potential can be flattened 
enough to support a stage of non-minimal inflation compati- 
ble with observations. The inflationary model was embedded 
in a moderate extension of MSSM augmented by three RH 
neutrino superfields and three other singlet superfields, which 
lead to a PQPT tied to renormalizable superpotential terms. 
The PQPT follows nMCI and resolves the strong CP and the 
/i problems of MSSM and also provides RH neutrinos with 
masses lower than about 10 12 GeV. The possible catastrophic 
production of domain walls can be eluded by the introduction 
of extra matter superfields which can be chosen so that the 
MSSM gauge coupling constant unification is not disturbed. 
For G masses larger than 8 TeV, observationally safe reheat- 
ing of the universe with T r h ~ 10 8 GeV can be accomplished 
by a three-body decay of the inflaton. The subsequent out-of- 



equilibrium decays of the produced RH sneutrinos can gen- 
erate the required by the observations BAU consistently with 
the present low energy neutrino data, provided that the Dirac 
neutrino masses are constrained in the range (1 — 10) GeV 
for all the light neutrino mass schemes. It is gratifying that 
the degeneracy of the masses of the RH (s)neutrinos required 
by the mechanism of nTL in our model is low enough com- 
pared with their decay widths, so that perturbative calculation 
remains safely valid. Finally, we briefly discussed scenaria in 
which the potential axino and saxion overproduction problems 
can be avoided. 
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